Exact three-body local correlations for excited states of the ID Bose gas 
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We derive an exact analytic expression for the three-body local correlations in the Lieb-Liniger 
model of ID Bose gas with contact repulsion. The local three-body correlations control the ther- 
malization and particle loss rates in the presence of terms which break integrability, as is realized 
in the case of ID ultracold bosons. Our result is valid not only at finite temperature but also for a 
large class of non-thermal excited states in the thermodynamic limit. We present finite temperature 
calculations in the presence of external harmonic confinement within local density approximation, 
and for a highly excited state that resembles an experimentally realized configuration. 
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When ultracold bosons are confined to move in only 
one dimension (ID), they provide a very clean realiza- 
tion [IH3] of a seminal exactly solvable model introduced 
by Lieb and Liniger (LL) pQ. Being an integrable model, 
it has very special dynamics showing almost no relaxation 
in experiments [3]. 

This fact stimulated lots of theoretical interest to un- 
derstand the thermalization of isolated ID systems and 
the role of integrability as well as its breaking in this pro- 
cess [5]-[9]. In particular, it has been shown [El [10], that 
virtual excitations of bosons to higher transverse modes 
of a confining potential result in a weak three-body local 
interaction that violates integrability of the many-body 
problem. Thus it is important to understand three-body 
local correlations in the absence of integrability breaking 
terms first. Such correlations have also been measured 
recently using analysis of particle losses [TTJ [T2] , density 
fluctuation statistics [13] , time-of- flight correlation statis- 
tics [13], and scanning electron microscopy [TS]. They 
provide a very sensitive test of coherence, and e.g., for 
Bose-Einstein condensates they increase by a factor of 
6 = 3! if the temperature is raised to be much larger than 
the condensation temperature [THJ E] . In spite of the LL 
model being integrable, analytical calculation of its corre- 
lation functions is notoriously hard [TS] . Two-body local 
correlations in equilibrium can be simply obtained using 
the Hellmann-Feynman theorem and exact thermody- 
namics, and show excellent agreement with experiments 
[T§] . Three body local correlations were analytically cal- 
culated only at zero temperature in a remarkable tour de 
force [20., as well as numerically in Ref. [21) . 

In this Letter, we exactly evaluate three-body local 
correlations in the thermodynamic limit for a large class 
of excited states which can be described by density ma- 
trices diagonal in the energy representation. In particu- 
lar, we apply our method at finite temperatures and for 
highly excited states similar to the ones created in exper- 
iments [3], and we take into account external harmonic 
confinement within local density approximation (LDA). 
We note that local two-body correlations in ID play the 
role of the "contact" introduced by S. Tan [521 123] ■ Simi- 
larly, three-body local correlations correspond to a three- 



body contact which is being actively explored [24 . Our 
exact results provide an important benchmark for such 
theories, as well as for numerical methods for simulating 
field theories in ID [55] . 

The model. — The LL model describes a system of 
identical bosons in ID interacting via a Dirac-delta po- 
tential. The Hamiltonian in second quantized formula- 
tion is given by 

H = [ dx^- (d x ^d x ip + cil>^iH>) , (1) 
Jo ^ m 

where c > in the repulsive regime we wish to study, 
and m is the atomic mass. The dimensionless coupling 
constant is given by 7 = c/n, where n — N/L is the 
density of the gas. We will express temperature T in 
dimensionless units r = T/Td, where Td = h 2 n 2 /(2mk&) 
is the quantum degeneracy temperature. 

The exact thermodynamics of the model can be ob- 
tained via Bethe Ansatz [H [TS] . Each eigenstate of the 
system with N particles on a ring of circumference L is 
characterized by a distinct set of quantum numbers {Ij} 
that are integers (half- integers) for ./V odd (even). The 
wave function can be expressed in terms of N quasimo- 
menta {pj} that satisfy a set of algebraic equations 

N 

LPj + Y,0(Pj-Pk)=2Kl j , (2) 
fe=i 

where 6(p) = 2 arctan(p/c). The wave function is iden- 
tically zero if any two of the {Ij} coincide, which is 
reminiscent of the Pauli principle for fermions. In the 
Tonks-Girardeau (TG) limit c — > 00, {Ij} correspond to 
the quantum numbers of occupied single-particle states 
of free fermions. 

In the thermodynamic limit, if one wants to consider a 
mixed state diagonal in the energy basis, this is achieved 
by introducing a filling fraction < // < 1 in the space 
of quantum numbers, which plays a role similar to the 
occupation number of free fermions. All results of the 
present Letter are valid for // which have a finite ther- 
modynamic limit at constant I/N; the limiting function 
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should be piecewise continuous and normalized. For cal- 
culations, it is more convenient to define a function f(p) 
in terms of the quasimomenta: denoting by p(p) the max- 
imal allowed density of quasimomenta in the vicinity of 
p, the quasimomenta density for a mixed state is given 

by f(p)p(p)- 

Since all quasimomenta are coupled to each other by 
Eq. pi), the density p(p) is not independent of f(p): it 
satisfies the integral equation and normalization condi- 
tion 



g,(Y>^) 



pip) = 



2tt 



dp' 
2^ 



f(p')<p(p-p')p(p'), (3a) 



n= dpf(p)p{p) . 



(3b) 



with the kernel ip(p) = 2c/ (p 2 + c 2 ). In thermal equi- 
librium, f(p) has to satisfy a set of nonlinear integral 
equations [26 ,[27], but our results will be valid for more 
general f{p). 

Local correlations. — The local fc-body correlation 
functions are defined as 



(4) 



The first two of them are relatively easy to calculate: 
g% = 1, while 172 in equilibrium is given by the Hellmann- 
Feynman theorem [19l [27] . 

Here we report the results for k = 2 and k — 3 for 
general f(p), which can be written in terms of functions 
hmip) { m — 1?2) satisfying the following integral equa- 
tions: 



h m (p)=p m + 



dp' 
2n 



f(p'M P ~p')h m (p'). (5) 



In the case of 52 the final formula is 



92(1, r) 



2j 2 f dp 
2tt 



f(p) [27Tp(p)p 2 -h 1 {p)p] , (6) 



which agrees with the result of the Hellmann-Feynman 
theorem for thermal equilibrium 27J, but is more general. 
Similarly, for k = 3 the final expression is given by 



7 

53(7, r) = 



dp 

2n 



(p 4 + cV)27r / 9(p)- 



(4p 2 + (1 + 2/ 7 )c 2 ) ph x (p) + 3p 2 h 2 (p) 



2 7 3 



dpf(p)p(p)p 



(7) 



In the case when f(p) is even, in equilibrium for example, 
the last term in Eq. (j\ is zero because the integrand is 
odd in p. Both Eq. (|6j) and Eq. Q are Galilean invariant 
expressions [57] . In the following we will first consider an 
equilibrium case and then will proceed to highly excited 
states. 




FIG. 1: Local three-body correlator 53(7, r) as a function 
of dimensionless coupling 7 = c/n for the uniform system 
at fixed dimensionless temperature r = T/To, where To = 
h 2 n 2 /(2mkB) is the quantum degeneracy temperature. The 
inset shows the large 7 asymptotic behavior on a log-log scale. 



In Fig. [T] we plot the result in thermal equilibrium 
for fixed r as a function of the coupling 7. In partic- 
ular, at zero temperature our result agrees with that of 
Ref. [50], up to the precision of the numerical evalua- 
tion of both expressions, 10~ 3 . The behavior of c/3 is 
qualitatively similar to that of g 2 analyzed in Ref. [19] 
and it distinguishes three different physical regimes: (a) 
7 > max(l, \/t), strong coupling (TG) regime, #3 <C 1; 
(b) t 2 < 7 < 1, quasicondensate regime, 53 ~ 1; (c) 
7 < min(r 2 , y/r), decoherent regime, 53 rs 6. In the in- 
set of Fig. [T]the large 7 asymptotics are plotted together 
with the analytic forms of Ref. [19]: 53 - 16tt 6 /(157 6 ) 
for r = and g 3 ~ 9r 3 /7 6 for 7 2 >• r ^ 1. 

Harmonic traps. — Next we turn to the experimen- 
tally more realistic case of atoms confined in a waveg- 
uide with a harmonic longitudinal potential. The ID 
regime is reached if /1, k^T <C huj±, where p is the 
chemical potential and uj±_ is the transverse oscillator 
frequency [28 . If the density profile in the trap varies 
smoothly, the correlations can be calculated by combin- 
ing our exact results with LDA [25]. The relevant prop- 
erties of the gas can be characterized by the LL coupling 
70 and the temperature parameter To at the center of 
the trap. In Fig. [2] we plot the three-body correlator 
53(70, To) at the trap center and the normalized average, 
53(70,70) = /dx(-0 t3 (x)V' 3 (a;)) / (Jdxn 3 (x)), against 
the dimensionless temperature To for different fixed val- 
ues of 70. Similarly to the results of Ref. [29] for the case 
of <?2, we find that unless the coupling 70 is very small, 
53(70, t ) « 53(70, t ) at any temperature. 

The curves in Fig. [2] are related to the observed change 
in time of the particle loss rate in Ref. [T2] • With in- 
creasing temperature, the three-body correlations grow 
according to our result, which leads to a higher proba- 
bility of inelastic three-particle processes which in turn 
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FIG. 2: Local three-body correlators 33(70, to) averaged over 
the trap (dots) and 33(70, to) in the center (solid lines), as 
a function of the dimensionless temperature to for fixed di- 
mensionless coupling 70 in the center of the trap. At high 
temperature, 33 approaches the value 6 = 3!, reflecting Bose 
statistics. 
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FIG. 3: Local three-body correlator <?3 X (7,pi) as a function of 
pi for different values of 7. The excited state is characterized 
by the inner "Fermi quasimomentum" p\. The horizontal 
lines correspond to 3f :E (7,oo) = [33(27) + 932(27)] /4. The 
inset illustrates a typical quasimomentum density for 7 = 1, 
pi = 0.5c. 



raises further the temperature of the gas. This positive 
feedback causes a non-trivial dependence of the particle 
loss on time, and a detailed analysis of heating mecha- 
nisms is needed to describe the time dependences of the 
loss rates. 

Three-body correlations for highly excited states. — 
Since Eq. ([8| is valid for general distributions f(p) [3D], 
we can use our results ([6|,([7]) in situations where the sys- 
tem is neither in equilibrium nor is in its ground state. 
We will illustrate this by considering a state which is mo- 
tivated by the experiment of Kinoshita et al. [3J , where 
each atom was put in a momentum superposition state, 
after which the two clouds performed many oscillations 
without observable thermalization. The state created in 
the experiment is not an eigenstate, and the harmonic 
trap might play an important role. However, let us con- 
sider here a simple "caricature" eigenstate which might 
capture the behavior of 33, and hence the role of inte- 
grability breaking, in these experiments. This state is 
characterized by an f(p) consisting of two disjoint rect- 
angular "Fermi steps" symmetric with respect to p = 
at zero temperature: f(p) = 9{p 2 — pf) — 9(p 2 — p 2 ,) 
with P2 > pi > (see inset of Fig. [3|. We are inter- 
ested in the dependence of <?§ x (7, pi) on the "inner Fermi 
quasimomentum" p\ . If p\ is fixed then the "outer Fermi 
quasimomentum" , p 2 , is determined from the normaliza- 
tion ( 3b I . The momentum kick in Ref. [3J corresponds 
to p\/c of order one. In Fig. [3] we plot g§ x (7,pi) for 
fixed values of 7 as a function of p\. We find that the 
correlations grow with the momentum of the kick and 
they can become greater than 1. For large p\, the quasi- 
momentum distributions of left and right goers become 
approximately independent of each other. However, to 
obtain the correct limit as p\ — > 00 one needs to take 
into account deviations of 0{2p\/c 1) in Eq. ([3| from 



7T. This results in gf{~f,oo) = [93(27) + 9 52 (2 7 )] /4; in 
particular, g^ x (j —> 0,oo) = 5/2. Similarly, g^("f, 00 ) — 
[ 52 (2 7 ) + 2]/2, and gf (7 -> 0, 00) = 3/2 [27]. 

Derivation of Eqs. ^ and Q . — In Ref. [3T] a novel 
method was proposed to calculate the gk correlators 
based on the observation that the LL model can be 
viewed as the combined non-relativistic, weak coupling 
limit of the sinh-Gordon model. The resulting formula 
reads as [3DJ [3T] 



9k 



s— k 



n^/fe)^^!,.-.^), (8) 



2tt 



where the form factors Fg k \pi, . . . ,p s ) are the infinite 
volume s-particle diagonal matrix elements of the oper- 
ator il)^ k ip k , which can be obtained from known sinh- 
Gordon form factors [32 34 . These series were investi- 
gated previously by truncating them after the first few 
terms [3T]. Here we resum these series to all orders ob- 
taining closed analytical expressions for the local corre- 
lations. 

It has been proven in Ref. 30 that 

1 y] <^Oi2MP23) • ■ • <p(p a -i,s) , (9a) 



F i 2) = ~3 X! ^O^Mto) • ■ • <P(P,-I,s) Pi, 



(9b) 



and based on evaluations performed in Mathematica for 
the first few we conjecture 

F ^ = ~5 H ¥>Ol2MP2 3 ) • ..¥>(P 8 -1, 8 )X 



2*1,* [Pl,.-(Pl2+P23- 



+ Pl- hs )] , (9c) 



4 



where Pij = pi — pj and J2p denotes a sum over all 
permutations of {pj}. Below we will illustrate how series 
^ can be analytically resummed for <?2, and details of 
similar calculations for g 3 are presented in EPAPS [27] . 

We will use abbreviations dp = dp/ (2?r)/(p) and ipy = 
ip(pij). Using the symmetries of the integrand in Eq. Q, 
we have for c 3 <?2(7, t)/(27 2 ) 



[dpiPl J dp 2 cp 12 + J dp 2 J dp 3 (f 12 (p 23 + 

Jdpipi J dpi (fii2 p 2 + Jdp 2 Jdp 3 (pi 2 (f 23 P3 4 



dpi p\ ftnpipx 



dpi Pi [hi(pi) -Pi] , (10) 



where the terms in the first square bracket coincide with 
the iterative solution of the integral equation ( 3a ) . Simi- 



larly, comparison of the terms in the second bracket with 
the iterative solution of Eq. ^ leads to h\{p\) —p\. Now 
the second terms in the parentheses cancel each other and 
we obtain Eq. ([6f. 

In summary, we derived an exact formula for the local 
three-body correlation in a paradigmatic system, the ID 
Lieb-Liniger Bose gas. Given that exact expressions for 
correlation functions are scarce even in integrable mod- 
els, we emphasize the analytic nature of our result. Our 
non-perturbative formula is valid at any temperature for 
arbitrary value of the coupling 7, and for a large class 
of excited states of the system. The result can open the 
window to an analytic treatment of integrability break- 
ing perturbations and thermalization in nearly integrable 
systems. 

We acknowledge funding from The Welch Foundation, 
Grant No. C-1739, from the Sloan Foundation and from 
the NSF Career Award No. DMR-1049082. 

Note added. — While this paper was under review, 
we learned about Ref. |35j were some of our results were 
rederived and generalized. 
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Supplementary Material for EPAPS 
Exact three-body local correlations for excited states of the ID Bose gas 

THERMAL EQUILIBRIUM 

In thermal equilibrium |26j . the filling fraction can be written as f(p) = (1 + e £ ( p ))~ 1 , where the pseudo-energy e(p) 
satisfies the nonlinear integral equation 

*' = -^ + 2S-/^*-"')^( 1 + e -"''). <") 
and the free energy is given by 



F = L Ln-k B T J ^ log(l + e- £ W) 



(12) 



Eq. (11) is coupled to the equations for the density 



p(p) = ^- + I % Hp') <p(p - p') p(p r ) , (13a) 



n = J dpf(p)p(p), (13b) 
At T — 0, Eqs. (fl~3|) and ( 11 ) decouple and 1/(1 + e E ^) becomes a Fermi step function f(p) = 8(p 2 — p F ) where pf 



is determined from J^p F dpp(p) = n. 



DIMENSIONLESS VARIABLES 



It is useful to introduce dimcnsionless quasimomenta q = p/c, and to change slightly the notation for the densities: 
p{cq) — > p{q), and similarly for f(q) and e{q). Then the equations for the density become 

Pit) = ^ + / % fW) <p(q - </) pW) , (14a) 

1 = / dqf(q)p(q), (14b) 



7 

where the kernel is f(q) — 2/(g 2 + 1). In thermal equilibrium, f(q) = (1 + e^ 9 -*) -1 , and the equation for e(q) is 

s(q) = -a+ q ^- j^viq-q') log (l + e- £(9,) ) , (15) 
The dimensionless chemical potential a = p,jk B T needs to be determined by the self-consistent solution of Eqs. 



(14 1,(15). Once e(q) is found, the free energy F is given by 

F = Lr^t fa- 7 /^log(l + e-^))) 
2m \ J 2ir J 



(16) 



The dimensionless form of the functions /i m are defined as 



/<„,(</> = '/" + / vWM<7-<7')M<7')- (17) 



HELLMANN-FEYNMAN THEOREM 



For completeness, in this section we derive a formula which was established earlier in Rcf. |31j based on the 
Hcllmann-Feynman theorem. This theorem states that 

§(£)-<**»>-»£(£), <«> 
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where the free energy can be calculated from Eq. (16). In dimensionless variables, 



52 ( 7l r)=rA L- 1 log(l + e- £ W) 



(19) 



Now we recast it in a form which does not exhibit explicit derivatives. We will need the derivative of Eq. ( 15 1 with 
respect to 7: 



de(g) 
d7 



2«7 2 7 [dq e'(q) 



(20) 



We start by writing 



«.<7.t)-™--t/£ 1*1 + + (21) 

where the prime denotes derivative with respect to 7. The key step is substituting for 1/2-k in the last integral the 
rest of Eq. (fl4a|) : 



1 + e 5 ^) J 1 + e £ ( ? ) 



dqp(q) 



1 + e 6 ^) 



dq 



p(q) f d ? 



1 + e £ (9) 



(22) 



In the last step we used the fact that the kernel is an even function. Now we use Eq. ( 20 1 to express the convolution 
term and continue with the equalities as 



dq e'(q) 



2ir 1 + e £ (9) 



dq p(q) 

a' 27 
— + — 

7 T 



- / dq 



I _|_ e s(q) 
A - P(9) -2 

dq iT^ q ' 



p(g) 

1 + e £ (9) 



e (g) + a 



(23) 



where the terms with e' dropped out and we used Eq. (14b). Plugging this result into Eq. (21) even the a' terms 
cancel and we finally arrive at 



1 + e £ (9) 

We can rewrite the second integral using partial integration: 

^) = 2 1 *jdq Y M W) {-r 



IglogCl + e-W). 



dq 1 de(g) 

2^ 9 1 + e £ («) dg 



Differentiating Eq. (15) with respect to g we obtain 

de(q) 2 7 2 , 
dq t 



dq 1 -^de(g) 
2^TTe^^ q - q) ^f' 



(24) 



(25) 



(26) 



which shows that de(q)/dq satisfies the same integral equation as 2j 2 /rhi(q) (see Eq. (17)). Thus expression (25) is 
clearly the special case of the general result 



ff2 ( 7 ,r) = 2 7 2 / dqf(q)p(q)q 2 -2j 2 



dq 
2n 



f(<l)hi(q)q, 



(27) 



for equilibrium when /(g) = 1/(1 + e 6 ^). 

In the original dimensionful variables one also has to take care of the coupling constanst dependence of the kernel. 
It is interesting to note that the analogous derivation for the T = case is more subtle because then one has to 
differentiate the boundary of integration explicitly. 
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SUMMING UP THE SERIES |g) FOR k 



Every permutation in the form factors ( 9c ) contains all the momenta and not only the first and the last few of 
them, which seems to render the procedure used in the main text for k — 2 unfeasible. However, this difficulty can 
be overcome by rewriting the form factors exploiting their special structure: 



F s 3) = \ <^(<7l2M<723) ■ ■ • <p(Qs-l,s) 



it 



s-1 

1\,s ^ qd,d+i 

d=l 
s-1 



V^d, 



- 1 



rf+lj 



\ X! i ^(912)^(923) • • • <p(q s -i, s ) (<?ts + <ii,s) - 2 9i,s X! ^(112) ■ ■ ■ <p(qd-i,d)<p(qd+i,d+2) ■ ■ ■ v>{q s -x, s )qd,d+i 



(28) 



where in the last term the kernel <p{qd,d+i) is missing. We used here the explicit form of the kernel and X]d=i Qd,d+i = 
gi iS . In the last term in the second line the chain of kernels is broken and the polynomial only depends on momenta 
at the ends of the resulting chains, which allows us to use the same technique as before. 



First term 



Let us start with expanding the polynomial part 



qf, s + =(qf 



9?) 



-{^ql+qi)q s + iqWs- 
(9s + q 2 s ) - (M 3 S + q s )qi + 3q 2 s qj . 



(29) 
(30) 



The second line of the expansion can be obtained from the first line by the transformation qi o q si that is, by a 
reflection of the indices. Since the kernel f{q) is an even function, the form factor is invariant under this transformation 
(actually, every single term in the permutation sum is invariant). We are thus allowed to focus on the terms in the 
first line, taking the second line into account by multiplying the result by 2 which cancels the overall 1/2. 

• Terms containing q\ only. 
The subseries corresponding to the terms (qf + qf) is 



d<?i {qi + q\) ( / dq 2 / dq 3 ip{q 12 )ip{q 23 ) + / dq 2 / dq 3 / dq 4 <f{qi 2 )(p{q 23 )ip{q 3 4) + 



dqx{qi + q\) 



2np{q 1 ) - 1 - / dq 2 <p{qi 2 



where we used Eq. ( 14a 



• Terms containing q\,q s - 
For the term — (4<7i + qi)q s we have the series 



(31) 



dq 1 {4q 1 +q 1 ) / dq 2 / dq 3 ip{qi 2 )tp{q 23 )q 3 + / dq 2 / dq 3 / dq4,ip{qi 2 )ip{q 23 )^{q 3 i)q4 + 



where in the last step we used Eq. (17 1. 
• Finally, in the case of 3qf q$ we are similarly led to 



dgi(4g? +qi) 



hl{qi) -Ql~ dg 2 (p{q\ 2 )q 2 



(32) 



dqi^qx [ dq 2 dq 3 tp{qi 2 )ip{q 23 )q 3 + / dq 2 / dq 3 / dq 4 f{qi 2 )ip{q 23 )f{q 3 i)q 4 + 



3 / dqxqf 



h 2 {qi) - q\ - / dq 2 (p{q 12 )q 2 



(33) 
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Combining Eqs. ((3TJ33I wc find that several terms cancel each other and we arc left with 

dq(q 4 + q 2 )2irp(q)- f ~dq (4q 3 + q)h x (q) + 3 / dqq 2 h 2 (q) 

dqi / dq 2 ip(q 12 ) ((qf - Aq\q 2 + Zq\ql) + (q 2 - qi q 2 )) . (34) 



(35) 



Using the symmetry of tp(qi — q 2 ) and then its explicit expression we can write the last term in the nice form 



7. I d<7i / dq 2 ip(q 12 )(q 4 2 + qf 2 ) = -~ / dqi I dq 2 2q 2 12l 



so the total contribution from the first term of the form factor is 

dq{q i + q 2 )2wp(q)- / dq{4q 3 +q)h 1 (q)+3 I dqq 2 h 2 (q) - I d Ql I ~dq 2 q\ 2 . 



(36) 



The polynomial part is 



Second term 



Qi,sQd,d+i = (liQd - qiid+i) + (qd+iq s - qdq s ) 



(37) 



Again, the second two terms give the same contribution as the first two due to the sum over d, so we focus on the 
first two terms and multiply the result by 2 in the end. Due to the missing kernel, the s-fold multiple integrals split 
into the product of a c£-fold and an (s — (i)-fold integral. 

• -2qiq d . 

Taking into account all the prefactors, the subseries corresponding to the first term in the first parenthesis is 



2 / dq 1 q 1 [ / dq 2 / dq 3 ip(q 23 )+ / dq 2 / dq 3 / dq i ip(q 23 )ip(q 3i ) + 



2 J dqi J dq 2 ip(qi 2 )qiq 2 (J dq 3 + J dq 3 J dq 4 



<P(<?34) 



2 J ^ J ^ 2 J ^ q3(P ^ qi2 ^ q23 ^ qi<l3 ^ Q4 + J ^ J ^5 V{<14:5) + ■ ■ ■ j +•■• , (38) 

where we reshuffled the series: the first line contains the d = 1 terms, the second line contains the d — 2 terms 
and so on. Let us add and subtract the term 2 J dqi J dq 2 q\ from the first line. Then all the parentheses become 
equal to I/7 (c.f. Eq. (14b)). This can be factored out leaving another infinite series 
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d-QlQl dg 2 ^(912)92+ / dq 2 / dq 3 ip(q 12 )tp(q 23 )q 3 



dqi 91/11(91) ■ 



Thus the contribution of the first term in Eq. (37) is 

2 
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dq q hi (q) + 2 / dqi \ dq 2 q 1 



(39) 



(40) 



• 2q 1 q d+1 . 

We can proceed in the same way for the second term in the first parenthesis of Eq. (37). 



2 / dqi qi ( / dq 2 q 2 / dq 3 <p{q 23 ) + J dq 2 q 2 J dq 3 J dq 4 tp(q 23 )<p{q 3 4:) + ... J + 

2 J dqi q x J dq 2 ip(q 12 ) (Jdq 3 q 3 + J dq 3 J dq 4 q 3 (p(q 34 ) + ■■■ j + 
2 / dq 1 q 1 / dq 2 / dq 3 tp(q 12 )ip(q 23 ) ( / dg 4 g 4 + / dq 4 q 4 / dq 5 tp(q i5 ) 



(41) 
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where we reshuffled the series similarly to the previous case. Let us add and subtract now the term 
2 J dqi /d<72<7i<22 ; so that all the parentheses above become equal to J dq 2irp{q)q. After factoring these out 
we are left with the same infinite series, so the contribution of the second term in Eq. (37) is 



dq 2irp(q)q 



2 / dqi / dq 2 qiq 2 ■ 



Thus the total contribution of the second line of Eq. ( 28 1 is 

2 

7 



dqqh 1 (q) + 2[ dq2np(q)q 



dqi / dg 2 (2<?i - 2q x q 2 ) . 



(42) 



(43) 



The final result: closed expression for 33(7, t) 



The final result is given by the sum of the partial results ( 36 ) and ( 43 1 . The last terms of these expressions exactly 
cancel each other because 

dqi J dq 2 q\ 2 = J dft J dq 2 {2q\ - 2q x q 2 ) , (44) 
so we arrive at the dimensionless form of Eq. ^ 

93(l, t) _ f , „ 2 ^_„^ f 3„ ( A„2 , 1 , ^ \ „ u f ~\ i o f ~A„„2 U f„\ , o ( f j„o_„/„N^ 2 (45) 



dq (q 4 + q 2 )2np{q) - / dq [ Aq z + 1 + - ) qh X (q) + 3 / dqq z h 2 (q) + 2[ / dq2TTp(q)q 



GALILEAN INVARIANCE 



In this section we show that our expressions (27), (45) are invariant under a Galilean boost, which provides a 
non-trivial consistency check of our results. 



The integral equations after a boost with momentum b read as 

dq' 



h m (q) 



2tt 



f^-b)ip(q-q')h m (q') 



Note that as a special case ho(q) — 2np(q). After a simultaneous shift in both variables q and q' we obtain 

dq' 



h m (q + b) = (q + b) m + 



2tt 



f(q'Mq-q')h m (q' + b). 



From the iterative solution of these equations it is easy to see that 

+ b ) = ho(q) , 

h 1 {q + b) = h 1 (q) + bh Q (q), 

h 2 (q + b) = h 2 (q) + 2bh 1 {q) + b 2 h (q) . 



(46) 



(47) 



(48a) 
(48b) 
(48c) 



Now we can calculate the boosted version of g 2 . From Eq. (27) 

dq 



92 
2 



7 



2tt 



f(q-b) 2np(q)q'-f(q-b)h 1 (q)q 



Shifting the integration variables and using Eqs. (48) we arrive at 

h 92 , f dq 



sy2 ry2 



2tt 



f(q) [h (q)(2bq + b 2 ) - hx{q) b - h (q)(q + b)b] , 



(49) 



(50) 



where we used ho(q) = 2np(q). Now due to J dq f(q)ho(q) q = j dqf{q)h\{q) the integral can be readily shown to 
vanish, implying the Galilean invariance of g 2 . 

The invariance of (73 can be shown along the same lines. 
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LOCAL CORRELATIONS IN THE "CARICATURE STATE 



In this section we discuss the derivation of the asymptotic values of the local correlations for an infinite momentum 
kick. Let us denote the various functions in this excited state by the superscript "ex" . The dimensionless integral 
equations can be written as 



<l2 



91 cV 
2^ 



'12 



91 



2tt 



+ )<p(q- q')K*(q) = 9™ + 



Q2 



91 



cV 
2^ 



In particular, /i§ x ((?) = 2np cx (q) with normalization 

I" 12 dq 



9i 



[<p(q - q') + (-!)">(<? + q')}h%(q'). (51) 



(52) 



We want to use as a reference state the ground state solution on the interval [— (q 2 — 9i)/2, (q 2 — <Zi)/2] = [— d, d], so 
we change our variables as 



, . qi + 92 _ , . 

q = k A = k + a . 



(53) 



This leads to 



h^(k + a) = h m (k) =(k + a) m + 



{k + a) m + 



dk' 
2^ 

dk' 
2tt 



[cp(k - k') + (-l) m tp(2a + k + k')} h m (k') = 



ip{k-k') + 



(-1)" 
2a 2 



k + k' 3(fc + fc') 2 -l 
a 4a 2 



h m (k') (54) 



and J_ a d §h (k) = 1/(27). 

Expanding the resolvent of the integral equation in a^ 1 we obtain the solution 



( — ~\\ m f d dk 

hm (k) = h b m (k)+ ( —^(l-0)- 1 o 



2a 2 



2tt 



1 - 



k + k' 3(fc + fc') 2 -l 
a 4a 2 



h h m (k') + 



4a 4 7o 



dfc 
2^ 



where /i^fc) is the solution of the equation obtained by keeping only non-negative powers of a: 

h h m (k) = (l-0)- 1 o(k + a) m . 



h h m (k)+0(a- 5 ), 
(55) 

(56) 



These are nothing else but the Galilean boosted functions discussed in the previous section (c.f. Eq. (48)): 

h%(k) = h§\k), h\(k) = hf\k) + ah ° ) (k), h 2 J {k) = h { 2 °\k) + 2ah ( °\k) + a 2 h ( °\k) 1 (57) 



where the superscript denotes the ground state solutions on the interval [—d, d]. The last term in Eq. (55 1 is the result 
of multiple applications of the operator (1 — <^) _1 , and we introduced the notations 



-d 



dk 
2^ 



h£\k) 



1 

To 



dfc, 
27 



hf\k)k = ei 



Using Eqs. (57) for h^ n (k) we find from Eq. (55) 
1 



ho{k) = h h (k) 
h^k) = h\{k) 

h 2 (k) = h h 2 (k) 



2 7o a 2 u w 2 7o a 3 1 w 8a 4 

^<\k) + -^hf\k)-\ 
2 7 o a 270 a z a A 



3e 2 



1 



dfc 
27 



+ ^)h( °\k) + -h 2 »'(k) 



h { °\k)k 2 = e 2 . 



» 62 -2 ei 



7o 75 
1 1 1 
870 4 7 2 



7o 



+ 0(a- 5 ) 



(58) 



(59a) 



4 0) W + ^r4 0) W 

870 



+ 0(a- 4 ), 



2 la y > 2 7o a 
,(o) 



-*i 0) (*) + ^ 



1 1 

870 4 7 2 



4 0) (fc) + ^4 0) w 
°7o 



(59b) 
(59c) 



Comparing the integal of ft (fc) and ho(k) we can relate 7 and 70: 



70 = 27+^+O(a" 3 ). 



(60) 
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Now we are in the position to calculate the local correlations. For g| x (7, a — > oo) we only need the leading corrections 
to ho(k) and h\(k) and we obtain 

= 2 £ 2 g [h?( q )q 2 firm = *fj^ + af hl (k)(k + a)] = ^ + 2 L+0 (a - 1 ) , (61) 

which gives #2 X (7, a — > oo) = g 2 (2'f)/2 + 1. For (/3 X ( 7 , a — > oo) we need all the corrections given above. The divergent 
positive powers of a cancel, and after a lengthy calculation we arrive at the remarkably simple result 

9 r(l,a -> oo) = i fl3 (2 7 ) + ^ 2 (2 7 ) . (62) 

It is interesting to note that as 7 — > 0, the asymptotic values are given by g§ x (7 —¥ 0, 00) — > 3/2 and g§ x ( 7 — > 0, 00) — > 
5/2. On the other hand, in the free boson case we recover <7! x (0, 00) = 2/2 + 1 = 2 and gfi®, 00) = 1/4-6 + 9/4-2 = 6. 



